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Abstract 

We construct the exact ground state for an antiferromagnetic spin-| 
model on the two-leg ladder as an optimum ground state. The ground 
state contains a discrete parameter a = ±1 and a continuous parame- 
ter a which controls z-axis anisotropy. For most values of a the global 
ground state is unique. It has vanishing sublattice magnetization and ex- 
ponentially decaying correlation functions. By using the transfer matrix 
technique, we calculate exactly the fluctuations of the magnetization, the 
nearest-neighbour correlation, and the longitudinal correlation length as 
functions of the parameters. 

Dedicated to Prof. H. Horner on the occasion of his 60th birthday. 

1 Introduction 

The investigation of quantum spin ladders is a very active field of condensed 
matter physics. Reports concerning experimental realizations of such systems 
are contributed frequently, see for instance A review can be found in 

H . From a theoretical point of view, spin ladders have been investigated by 
using exact diagonalization W, density matrix renormalization group (DMRG) 
methods || ^ bosonization [|10|, and various other techniques, see e.g. |ll[] jl3| 
and references therein. 

Recently, Kolezhuk and Mikeska presented a set of matrix product 
ground states (MPG) for special isotropic spin-| ladder models. To each rung of 
the ladder, a matrix is assigned, which has local spin states for the corresponding 
rung as its elements. The global state is given by the product of these matrices, 
in which the matrix elements are multiplied via the tensorial product in spin 



space. If the model parameters obey the conditions given in 14 , the resulting 
global state is a so called optimum ground state, i.e. it is not only a ground state 
of the global Hamiltonian, but also of every local interaction operator. 

In the present work we construct a one-parametric set of optimum ground 
states for antiferromagnetic spin-| ladders with two legs. The local interaction 
is identical to the one on the hexagonal lattice in our previous paper jL5|. In 
contrast to Jl4| , the global ground state is given in terms of a vertex state model, 
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which is a generalization of the MPG approach. This allows to specify the local 
spin states for every lattice site instead for complete rungs as required by the 
MPG approach. Properties of the ground state are calculated by using the 
transfer matrix technique, which is explained in Appendix A. 



2 The model 



Consider a ladder of length L with two legs and periodic boundary conditions. 

3 
" 2 ' 



Each lattice site is occupied by a spin-|. The global Hanhltonian 



H 



L 

Ei 

i=l 



hi,i+l + hi 



(1) 



contains only nearest neighbour interactions. If the ladder is visualized hori- 
zontally, the index i (i 1 ) denotes the upper (lower) spin on rung number i. All 
local interactions are equal, they only act on different spin pairs, i.e. the system 
is completely homogeneous. 

The local interaction operator is the same as in the spin-| model on the 
hexagonal lattice presented in our previous work jlq ]. It is given in terms of 
projectors onto its eigenstates: 



where 
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A3, A 2 X±2, X 2 are positive real numbers, a is real, and a = ±1. The canonical 
spin- 1 basis states are denoted as 



S*|3) = 
S z \\) = 



113) 



S*|3) = 
S Z \T) = 

Note that (||) has the following properties: 



-113) 
-111)- 



(4) 



1. It has rotational symmetry in the xy-plane of spin space, i.e. it commutes 
with the local magnetization operator Sf + Sj . 

2. It is parity invariant, i.e. it commutes with the operator Py, which inter- 
changes the spins at lattice sites i and j. 



2 



3. It has spin-flip symmetry, i.e. it is invariant under the transformation 
S z -> -S z . 

4. Its lowest eigenvalue is zero, i.e. is positive semi-definite. 

The Hamiltonian contains 5 continuous parameters, namely 4 A-parameters plus 
a. This includes a trivial scale. The two-spin states (||) are the excited local 
eigenstates of hu, the remaining 9 eigenstates are local ground states, i.e. the 
corresponding eigenvalue is zero. 

Since the global Hamiltonian (]l|) is a sum of positive semi-definite operators, 
zero is also a lower bound of the global ground state energy Eg. In the next 
section a global eigenstate corresponding to eigenvalue zero is constructed, which 
must therefore be the global ground state. 

For o = —y/3 and a = — 1, the A-parameters can be chosen so that (p|) are 
eigenstates of (Si + Sj) 2 . In this isotropic case, the local interaction (g) has 
complete SO (3) symmetry in spin space and can be written as 

^ = ^-S, + ^(S,.S,) 2 + ^(S i .S,) 3 + ^ ! (5) 
which simply projects onto all states with (Sj + Sj) 2 = 3(3 + 1). 



3 The global ground state 

To each lattice site we assign a set of vertices with binary bond variables which 
are denoted as arrows. The values of these vertices are single-spin states at the 
corresponding lattice site. On the upper row 



<T> , x > T < 

V : a|3> A 
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On the lower row 
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a|3> 



\/ 
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The parameters a and cr are the same as in (Jf^j. 

In order to construct the global ground state from these local spin states, it 
is necessary to define the concatenation of the vertices (^) and (0) , respectively. 
This is similar to classical vertex models of statistical physics, but the generic 
product of numbers is replaced by the tensorial product in spin space: 



A V 



+ 




(8) 



i.e. the interior bond between the two vertices is summed out. Concatenations 
between adjacent vertices on the lower row and between upper and lower vertex 
on the same rung are completely analogous. The resulting cluster represents a 
state in the Hilbert space of two neighbouring spins. 

A global state \^q) can be constructed by successively attaching vertices to 
the cluster until the ladder has the desired length. This process is associative, 
i.e. the order in which the vertices are concatenated does not matter. As periodic 
boundary conditions are imposed, the 'free' bonds at the first and last rung of 
the ladder have to be summed out. Because of the similarity to classical vertex 
models, such a global state is called a vertex state model |l5|, [l7]]. These are 
direct generalizations of the well-known matrix product ground states for spin 
chains |16, If |(pf lM2M3 ) denotes the vertex at site i with arrow variables /ii, 
fj,2, and /i3 (cf. definitions (||) and (Q)), a formal expression for the ground state 
is given by 



l*o> 



e n w 1 ' 



(9) 
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The sum is over all arrow configurations on the bonds. 

In order to show that the resulting global state \^>q) is the ground state of 
the Hamiltonian (|l|), we collect all two-spin state which are generated by all 
possible concatenations of two vertices^]: 

|31)+ct|13) |3T)+a|l3) 

|ll)+o|3T) |TT)+a|31) 

|ll) + o|T3) |TT) + a|13) (10) 

|1T) + era 2 1 33) |Tl) + era 2 1 33) 

|lT) +C r|Tl). 

It is easy to check that each of these 9 two-spin states is perpendicular to all 
local excited states (|3|), i .e. they are annihilated by the local interaction operator 
(|[). In other words, (|Tc|) are the local ground states of hij. 

From the construction of the vertex state model it is clear that any projection 
of l^o) onto the Hilbert space of two neighbouring spins is a linear combination 
of the two-spin states ([To|). Therefore |^o) is annihilated by all local interaction 
operators, 

M*o} = o (ii) 

for all nearest neighbours i and j. Hence 

H|*o)=0. (12) 

As explained above, zero is a lower bound for the ground state energy Eq, 
therefore |^o) must be the ground state of H. This type of global ground 
state, which simultaneously minimizes all local interaction operators, is called 
an optimum ground state flTBt-fTijl, since E takes the lowest possible value. 



4 Properties of the ground state 

The first interesting expectation value is the magnetization of a single spin. As 
can be seen from the vertices @ and (Q), a spin flip S z — > — S z is equivalent 
to a flip of all arrows on the bonds, which leaves the global state invariant. 
Therefore the single-spin magnetization (Sf)\s must vanish. In this sense the 
global ground state is antiferromagnetic. 

For a — > oo, the structure of the global ground state becomes very simple. 
In this limit, the ground state is completely dominated by the vertices 
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(13) 



1 Common prefactors have been omitted. 
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Figure 1: Fluctuations of the pair magnetization as a function of the parameter 
a 



Clearly, there are only two different possibilities to concatenate these four ver- 
tices to form a global state. These two possibilities correspond to the two 
different Neel states on the ladder. Thus in the limit a — > oo, the global ground 
state is simply given by^| 



|*o) = |Neeli) +a L |Neel 2 ) 



(14) 



The transfer matrix technique can be used to calculate expectation values for 
arbitrary values of the parameters a and a. Appendix A explains how to apply 
this method to the two-leg ladder. For the present model, the eigensystem of 
the transfer matrix can be obtained exactly, but has a very complicated form. 
For this reason, explicit formulae for the expectation values have been omitted. 
Instead, the results for the thermodynamic limit are plotted as a function of the 
parameter a. The dependence on a drops out in all expectation values, which 
are calculated in this work. 

Consider the combined magnetization operator of an elementary cell at rung 
number i, 



Sf 



qZ 



(15) 



where the index iji' again denotes the upper/lower spin on the rung, respec- 
tively. Of cou rse, (Sf) oo vanishes. Its fluctuations can be calculated by using 
equation (A. 6). The result is shown in figure [l| The maximum at a rs 1.4 can 
be understood as follows: It is clear from the vertices @ and (0), that for a = 
the system contains only the single-spin states |1) and |1). With increasing 
values of a, |3) and |3) states are mixed in, so the fluctuations become larger. 
Finally, for a — ► oo, the global ground state is simply a superposition of both 
possible Neel states, so the system is frozen. Therefore, the fluctuations must 
vanish in this limit. 

Correlations between two different rungs can be obtained from equation (|A.7|) . 
The nearest neighbour correlation, (Sf 5 f f +1 ) o, is plotted as a function of a in 
figure |j. The shape is very similar to the one in figure [j]. Note that the correla- 

2 The global prefactor a 2L has been dropped. 
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tion is antiferromagnetic for all values of a. The vanishing in the limit a — > oo 
becomes clear from 

{SfSf + i)oo = wf^f+l/oo + (S^Sf +1 ,)oc 

The first two terms approach — | in the Neel limit a — ► oo, but the last two 
terms converge to +|, since the spin operators act on the same sublattice. 
The minimum at a « 1.5 is very close to the maximum of the fluctuations (cf. 
figure 0). 

As a function of the distance r, the correlation between Sf and S* decays 
exponentially, i.e. (S*S*) oo oc exp(— r/£j). Figure || shows the a-dependence of 
the corresponding inverse correlation length O^ 1 . Note the divergence at a = 1. 
At this point, the weight of all non-vanishing vertices of the corresponding 
classical vertex model (cf. Appendix A) is 1. This corresponds to an infinite 
temperature, if the vertex weight is interpreted as a 'Boltzmann weight'. The 
correlation length remains finite for all values of a, so the system is never critical. 
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For a ^ 0, the constructed optimum ground state |^o) on the ladder is the 
only ground state of the global Hamiltonian ([!]). The proof can be carried out 
by induction according to system size and is omitted here. In the special case 
a = 0, the global ground state degeneracy grows exponentially with system size. 



5 Summary and outlook 

We have constructed an optimum ground state for a 5-parametric spin-| model 
on the two-leg ladder. The Hamiltonian is completely homogeneous and contains 
only nearest neighbour couplings. The local interaction operator exhibits S z - 
conservation and parity invariance. No external magnetic field is applied. For 
special values of the parameters, the system has complete SO(3) symmetry. 

Optimum ground states are not only ground states of the global Hamiltonian, 
but also of every local interaction operator. For the present model, the ground 
state is given in terms of a vertex state model. A set of vertices, which have 
single-spin states as their values, is assigned to each lattice site. Concatenating 
these vertices on the two-leg ladder yields the global ground state. Vertex state 
models are straightforward generalizations of the well-known matrix product 
ground states for spin-chains. The ground state contains a continuous parameter 
a and a discrete parameter a = ±1. 

Properties of the ground state have been calculated by using the transfer 
matrix technique. The sublattice magnetization and the total magnetization 
vanish, so the ground state is antiferromagnetic. Two-point correlations along 
the ladder decay exponentially. The fluctuations of the sublattice magnetization, 
the nearest neighbour correlation, and the longitudinal correlation length have 
been determined as a function of the parameter a. As the correlation length is 
always finite, the system is never critical. Except for a = 0, the global ground 
state is unique. 

For other spin-| models on the two-leg ladder, e.g. with generic Heisenberg 
interaction, the constructed ground state could serve as a variational ground 
state. In this scenario, the anisotropy parameter a plays the role of a vari- 
ational parameter. Optimum ground states can also be constructed for more 
sophisticated ladder models, e.g. with next-to-nearest neighbour interaction and 
inhomogeneous couplings. 



Appendix A: The transfer matrix technique 

Consider a local observable A4, which acts on the spins on rung number i. Its 
expectation value in the ground state is defined as 

The denominator can be interpreted as two identical vertex state models on top 
of each other. The upper (lower) one represents the bra- (ket-) vector. Since the 
vertices at a given lattice site generate only local spin states, the inner product 
can be taken separately at each lattice site. The result is a classical vertex 
model, i.e. the 'weights' of the vertices are numbers, but each vertex emanates 
two bonds in each direction. One set of bonds carries the arrow variables of the 
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bra-vector, the other one carries those of the ket-vector. Therefore the inner 
product (v&ol^o) can be constructed from the L-fold product of the transfer 
matrix 



T, 



(^3,^4,1/3,1/4) 



0*1,1*2, "\,V 3 ) 
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(A.2) 



It has 4 incoming (^1,^2,^1,^2) and 4 outgoing indices (/13, /14, 1/3, 1/4). Each 
index represents a binary bond variable, so T is a 16 x 16-matrix. As we are 
dealing with periodic boundary conditions, the trace must be used to obtain the 
inner product 



(l-o |*o) =trT z 



(A.3) 



The numerator of (A.l) can also be interpreted as a classical vertex model. 
It differs from the one for the denominator only at rung number i, where the 
observable A; t acts non-trivially. Therefore the numerator also has a matrix 
product representation, but the i-th factor is the associated transfer matrix 
T{A) instead of T: 



(*o|^|*o)=tr [T^ 1 T{A)T L - 1 ] = tr [T^T^ 1 ] 



(A.4) 



T{A) is obtained from definition (A.2) by inserting Aj b etween the bra- and the 
ket-vector on the r.h.s. of the equation. Note that ( A.4 ) is independent of i, as 
l^o) has perfect translational invariance. Since T is symmetric, its eigenvalues 
Xk and its normalized eigenvectors \u k ) can be used to calculate the L-fold 
matrix products. This yields 



E k x L k 



(A.5) 



for the expectation value of Aj. In the thermodynamic limit L — » 00, which 
is the most interesting case, only the largest eigenvalue Xmax survives. This 
simplifies (A.5) to 



(A)oo — 



Xn 



(A.6) 



The above consideration can be easily extended to two-point correlation 
functions. In the thermodynamic limit the corresponding formula is 



{A 1 B r ), 



O 

A max 



\T(A)\u k )(u k \T(B)\u n 



Xk 



(A.7) 
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If the following mapping of matrix indices to arrow configurations is used, 



the transfer matrix (A. 2) takes a block-diagonal form 
index 



9 

10 



11 

12 
13 
11 



15 



upper bra 



upper ket 



lower bra 

A*2/M4 



lower ket 



(A.8) 
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The first block is the 6 x 6-matrix 
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(A.9) 



This submatrix contains the leading eigenvalue of the whole transfer matrix for 
all values of the parameters. 

The 4 x 4-blocks for indices 7 — > 10 and 11 — > 14 are identical. Both are 
given by 



Tr- 



io 



Ti 



1->14 



/ 2cr a(l + a 2 ) 
cr(l + a 2 ) 2cr 

1 1 
V 1 1 



1 
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2(7 

(7(l + a 2 ) 
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(7(l + a 2 ) 

2(7 



(A.10) 



The remaining 2-dimcnsional subspace of indices 15 and 16 is already diagonal: 
Xi5 = Xi6 = 2 (A.ll) 



Since both submatrices (A.9) and (A. 10) are symmetric, an orthogonal eigen- 
basis of the complete transfer matrix exists. However, exact expressions for the 
eigensystem are very complicated, so they have been omitted. 
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